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I Abstract. 

^ ■ We present ^^CO J=l-0 and J=2-l Swedish-ESO Submilhmetre Telescope (SEST) observat ions of the barred 

spiral galaxy M83 (NGC5236) . The size of the CO maps is 10' x 10' and they cover the entire optical disk. The 
grid spacing is 11" for CO(J=1-0) and 11" or 7" for CO(J=2-l) depending on the position in the galaxy. In 
^ ' total we have obtained spectra in 1900 and 2574 positions in the CO( J=l-0) and C0(J=2-1) lines, respectively. 

The CO emission is strongly peaked toward the nucleus, which breaks up into two separate components in the 
CO( J=2-l) data due to the higher spatial resolution. Emission from the bar is strong, in particular on the leading 
edges of the bar. The molecular gas arms are clearly resolved and can be traced for more than 360°. Emission in 
the inter-arm regions is detected. The average CO ( J=2-l)/( J=l-0) line ratio is 0.77. The ratio is lower than 
' this on the spiral arms and higher in the inter-arm regions. The arms show regularly spaced concentrations of 

O molecular gas. Giant Molecular Associations (GMA's), whose masses are of the order lO'^ Mq. The total molecular 



OO 

o 



gas mass is estimated to be 3.9x10^ M©. This mass is comparable to the total HI mass, but H2 dominates in the 
optical disk. In the disk, H2 and HI show very similar distributions, including small scale clumping. We compare 
the molecular gas distribution with those of other star formation tracers, such as B and Ha images. 

b , 

^ ^ Key words. Galaxies: individual: (M83; NGC5236) - Galaxies: spiral - Galaxies: structure - Galaxies: ISM - Radio 

C/3 . lines: galaxies - Galaxies: abundances 



1. Introduction role in understanding the structure and evolution of galax- 
j; ics. Unfortunately, direct studies (in emission) of the cold 
^ . In disk galaxies the mass of the molecular (H2) gas is usu- (^^^q k) H2 gas in molecular clouds are very difficult, since 
- - ■ ally similar to, or a bit less than, that of the atomic (HI) ^^^^^^ of ^^^^ ^^.^^ ^^^-^^^ rotational level is of the or- 
gas ( |Young fc Scoville| | l99lMCasoli et al. || l998D . However, ^qqj^ (corresponding to a wavelength of 28 ^^m), 
their spatial distributions are completely different, both therefore much larger than the average kinetic tem- 
in the radial and in the vertical directions. The molecular peratures in the clouds. Carbon monoxide (CO) is the 
gas IS concentrated to a relatively thm disk, which often ^^^^^^ ^^^^ common molecule with a relative abundance 
has a radial extent similar to that of the optical disk and ^-^^ ^^^^^^^ jj^ of ^^o^^ rphe heavy, asymmet- 
a vertical scale height which is smaller than that of the ^.^^^1 molecule has rotational levels which are easily 
young, massive stars. The radius of the HI disk is usually ^^^^^^^ ^y collisions with H2, and it radiates at mm wave- 
much larger than that of the optical disk (up to 5 times the jg^gf where observations can be carried out using ra- 
Holmberg radius), and the atomic gas is also less concen- telescopes. However, the conversion factor which is 
trated to the disk in the vertical direction ( |Knapp| | l987D . ^g^j ^^^-^^ ^.^j^^^^ density from the velocity- 
If the mass surface density in a region exceeds a cer- integrated CO intensity, Xco, is not well-known, and this 
tain critical dens ity, massive star formation is triggered is a source of major uncertainty ijYoung fc Scovillelll99ll: 
l)Kennicuttlll989^ . So, naturally, the study of interstellar |Combeslll99lh . 
gas, and especially the molecular gas phase, plays a key 

CO emission at mm wavelengths has been detected in 

Send offprint requests to: Andreas Andersson Lundgren (an- many galaxies ranging in distance from the nearby irregu- 

dreas@astro.su.se) lar galaxies The Large and The Small Magellanic Clouds 
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to quasars at a redshift of more than z~5, and in es- 
sentially all galaxy types. A large survey of 300 galaxies 
(o f which 236 were de tected) of different types was done 
bv I Young et al using the FCRAO 14 m telescope. 

The analysis of the data is still ongoing, but some in- 
teresting results have been published, e.g., the conversion 
factor (Xco) app ears relatively constant , at least within 
each galaxy type llYoung fc R,owndll200ll) . and the global 
star formation effi ciencY decreases with increasing galaxy 
size In this survey a number of galaxies 

were also mapped and radial distributions of the molec- 
ular gas can be examined. Recently, in a more detailed 
study, 28 nearby gal axies were m apped with the 45 m 
Nobe yama telescope l|Kuno et al ]|2000; N ishivama et al.l 
l200lll . Among other things, this showed that barred and 
non-barred galaxies differ in both the distribution and the 
kinematics of the molecul ar gas. Another large survey is 
the BIMA SONG project llRegan et al.ll200ll:lHelfer et all 
12003 ') . in which 44 nearby {vq < 2000 km s~^) galaxies 
have been observed with the BIMA array and the NRAO 
12m telescope. Early results show that many galaxies have 
multiple peaks at the nucleus, that the CO distribution 
usually cannot be fitted with a single exponential disk 
model, and that the ratio of the scale length in the CO 
disk and the scale length in the stellar disk is on average 
1 (but the scatter around this average is large). 

Detailed morphological studies of selected galaxies are 
more rare, partly because there are so few nearby galaxies 
for which the limited angular resolutio n of radio telescopes 
gives an adequate spatial resolution. iNakai et al ] il994) 
outlined a number of important observational conditions 
that should be fulfilled in order to be able to extract use- 
ful information on the properties of the molecular gas in 
different regions of a galaxy: (1) a spatial resolution bet- 
ter than one kpc, which is the typical width of molecular 
spiral arms and the spacing between them, (2) high sensi- 
tivity in order to detect the weak inter-arm emission, (3) 
mapping of a large field in order to get statistically mean- 
ingful results, and (4) filled aperture telescope data, since 
it is important to detect the diffuse emission in the inter- 
arm regions. To this list we add a fifth condition, (5) the 
observation should be carried out in (at least) two lines 
per molecular species (and preferably also partially ob- 
served in a few isotopomers) in order to gain insight into 
the density, temperature, and chemistry of the molecular 
gas. 

Single-dish measurements are therefore essential to our 
determination of the molecular gas distribution, its rela- 
tion to star formation activity, its location with respect 
to dynamical resonances, and the total reservoir of gas 
capable of forming future generations of stars. The most 
detailed single-dish maps of molecular gas in spiral galax- 
ies are the mappings of M 51, us ing the IRAM 30m tele- 
scope arcia- Burillo et al."l993') and the Nobeyama 45m 
telescope [Nakai ct a l. 1994), and tw o projects currently 
underway to map our neighbor M31 l|Loinard et alJll999l: 
iNieten et al.l200ffl . In M 51 the molecular spiral arms were 
resolved, and streaming motions were detected. Also, the 



CO emission was found to peak on the inside of the spi- 
ral arms, giving support to a scenario where molecular 
gas agglomerates and forms new stars, which in their turn 
produce Ha emission and other star formation tracers, 
close to the spiral arms. Support was also found for the 
proposal that gravitational instability in the disk is the 
main formation mechanism of GMAs (Galactic Molecular 
Associations of mass ~10'' Mq). Other large spiral galaxies 
that have been completel y surveyed in CO with a single- 
dish telescope are IC 342 jcrosthwaite et al.ll200lh . M 10 1 
llKennev et aLlll99lh. NGC 253 I Hoiighton et alJ Il997l) . 
and NGC 6946 l)Tacconi fc Young|ll989(l . Some large scale 
CO mapping proj ects, such as the OVRO ma ps of M51 
l(Aalto et alJl999l) and M 83 l|R,and et alJl999l) . have been 
done using mm-wave interferometers. Although this gives 
a superior angular resolution and important unique in- 
formation, a significant part of the emission is resolved 
out (especially the diffuse component in the inter-arm re- 
gions). 

We have mapped the CO(J=1-0) and C0(J=2~1) 
emission in the barred, grand-design, spiral galaxy M 83 
located at a distance of 4.5 Mpc. The observations were 
done with the 15 m Swedish-ESO Submillimetre Telescope 
(SEST) providing beam widths of 45" (1-0) and 23" (2-1), 
which corresponds to 980 pc and 500 pc, respectively, at 
the adopted distance. M 83 is among the brightest galax- 
ies in terms of CO emission, and it is also one of the most 
nearby barred spirals. Its low inclination (Ri24°) and prox- 
imity make it a perfect candidate to investigate the corre- 
lation between various tracers of stellar activity. 

In this paper we present the CO(J=1-0) and 
C0(J=2-1) observational data. We show the distribution 
of the molecular g inferred from the CO data, and 
compare it to the distributions of optical light, obtained 
in different filters, and the HI gas. We examine the CO 
( J = 2 — l)/( J = 1 — 0) line ratio and its variation over 
the disk. In two forth-coming papers we will present and 
discuss the kinematics of the molecular gas, and discuss 
the relation between the molecular gas and star formation. 

2. The galaxy M83 

M 83 is a relatively nearby, barred, grand-design, spiral 
galaxy (see Table QJ. It is viewed almost face-on and the 
bar is aligned with the major axis. Color images show 
clumpy, well-defined spiral arms, evidently rich in young 
blue stars. It is also fairly symmetrical and has no nearby 
massive optical companions and no immediate evidence of 
interaction or outflows. 

Due to its favorable properties i t is one of 
the most well-observe d galaxies: HI ("Rogstad et al.| 
llluchtmeier fc BohnenstengeL1981;.Tilanus fc Allen 



1974|:fHuchtmeier fc Bohnenstengell981; Tilanus fc AUer 
1^2^ ). radio continuum l(OndrecherJll985l lNjinmgg^^ 



a 



1993() , CO (see Tabled, I R continu um (Ad amson et al 
<198';ilFitt etal...l992; Elmepeen et al. 1998). optical con 



tinuum ll Jensen et alJll98lh. Balmer lines jComtd Il98lt 
Ide Vaucouleurs et al."l98r{': 'Tilanus fc Allen"l993^. and in 
X-ray emission (^Immler et al...l999.: .Soria fc Wu..2002i) . 



A. A. Lundgren et al.: Molecular gas in the galaxy M83 



3 



Table 2. Previous CO observations of M83 divided into two groups: single-dish and interferometer observations. If 
no location is given, the map is centered on the nucleus and oriented along the major axis. 



Reference 


Transition 


Beam and Telescope 


Map description 


Combes et al. CigTS) 


CO 


[1-0) 


a a" 

64 


NRAO 11m 


1 n nnsitinns fi' x 


Lord f 19871 


CO 


'1-0) 


45 


ICRAO 14m 


21 positions, 5'x5' 


Wiklind et a] ('19901 


CO 


^1-0) 


45 


bbal 15m 


_L \J JJVJOl IflVJllO, -1— J Vj^ U.Cl*-l-i- Ct±l Li ^ O /\ ^ 




CO 


1 2-1) 


23 


bbal 15m 




Handa et al (19901 


CO 


Jl-O) 


16 


NRO 45m 




Wall (19911 


CO 


'2-1) 


22 


JCMi 15m 




CO 


'3-2) 


22 


CSO 10m 




Petitnas Wilson (19981 


CO 


1 3-2) 


1 /l" 

14 


JOMl 15m 


81 positions, 0'7x0!7 


CO 


^4-3) 


1 o" 

12 


JOMl 15m 


46 positions, 0'7x0!5 


Israel fc Baas (20011 


uvj 


1 \ 
^Z-i ) 


zi 


jLyivii iom 


49 positions, l.'2x2!0 






'9 o^ 

^o-z ) 


14" 


JCMT 15m 


55 positions, l.'2xl!7 






'A Q\ 

^4-o ) 


11" 


JCMT 15m 


20 positions, 0.'5x0!5 


Crosthwaite et al. (20021 


CO 


Jl-0) 


55" 


NRAO 12m 


On-the-fly, lO'xlO' 




en 


'2-1) 


28" 


NRAO 12m 


On-the-fly, 8' x 8' 


Thuma et al. (in prep.) 


CO 


;3-2) 


23" 


SMT 10m 


374 positions, 5^5x2.'5 




CO 


;4-3) 


17" 


SMT 10m 


21 positions 


This paper 


CO 


;i-o) 


45" 


SEST 15m 


1900 positions, lO'xlO' 




CO 


;2-i) 


23" 


SEST 15m 


2574 positions, lO'xlO' 


Lord fc Kennev (19911 


CO 


;i-o) 


12'.'4x5'.'4, OVRO 


Inner eastern arm, l.'5xlf5 


Kennev & Lord (19911 


CO 


;i-o) 


8"9x5'.'8, OVRO 


Western bar end, lf5xl!5 


Handa et al. ( 19941 


CO 


;i-o) 


12" 


x6", NMA 


Nucleus, l!lxi:i 


Rand et al. ( 19991 


CO 


;i-o) 


6'.'5x3'.'5, OVRO 


Inner eastern arm, 2.'9xl^l 



Table 1. General Parameters of M83. 



Morphological Type" 


SAB(s)c 


IR center:* 




R.A. (J2000) 


13''37™00!8 


Decl. (J2000) 


-29°51'56" 


LSR systemic velocity (opt)'^ 


506 km s"^ 


Distance'' 


4.5 Mpc 


Position angle'^ 


45"= 


Inclination-^ 


24'' 


Holmberg diameter (Do)'' 


14f6 


Mm" 


7.7x10^ Mq 




3.9x10'' Mo 



^^VmjTOideu^^^l^ ^^7^; ^ofu^^^Wakanmtsu 

(19g4); '= Comte (19811 : Thim et al. (20031 : 

Huchtm eier fc Bohnenstengel (19811 : ^ Talbot et al. (19791 : ^ 
this paper 

The distance to M 83 is a long-standing issue. 
ISandage fc TammannI l)l974|l used a relation between the 
luminosity class and the angular size of the three largest 
Hll regions to deduce a distance of 8.9 Mpc, while 
Ide VaucouleursI |.197{/) used a number of methods to de- 
rive a distance of 3.7 Mp c. Recen t ly, Ce pheids have been 
observed in M83 by T him et all l|2003^ using the VLT. 
Twelve Cepheid candidates were observed and the dis- 
tance was estimated to be 4.5±0.3 Mpc. This distance 
is also supported by studies that seem to indicate that 
M 83 interacted w ith NGC 525 3 some 1-2 x 1 0^ years ago 
l|Kobulnickv fc SkiUman 1995t ICalzetti et al.l[l999.1 . This 
small, metal-poor dwarf galaxy about 2° SE of M 83 has 
roughly the same systemic velocity, and shows a highly dis- 
turbed HI velocity field and intense star formation. There 



has been a number of distance measurements to NGC 5253 
l|Saha et al . 1995: Parodi et al. 2000; Gibso n et al...2000|l . 
and they all favor a distance of about 4 Mpc. 

3. Observations and data reduction 

3.1. Observations 

The CO(J=l-0andJ = 2-l) observations of M 83 were 
done using the 15 m SEST^ during two epochs, 1989-1994 
and 1997-2001 [a d escription of the telescope is given by 
iBooth et all l(l989ll ]. Between 1988 and 1994 we observed 
CO(J=1-0) in the inner disk of M 8 3. A segm ent of these 
observations were presented in Wi klind et al . (1990). In 
1998 and 1999 we completed this map by observing the 
outer disk. During the latter observations it was possible 
to do simultaneous observations in two frequency bands. 
We used this option to observe the C0(J=2-1) line in 
parallel with the CO(J=1-0) line. The map spacing (11") 
and integration times were set by the requirements on the 
C0( J=l-0) data. In 2000 and 2001 our primary task was 
to fill in the central area of the C0( J=2-l) map. For these 
observations we chose a grid spacing of 7". We used the 
other receiver to observe the ^'^C0( J=l-0) line, but these 
data will be presented elsewhere. 

The receivers were centered on the CO(J=1-0) and 
C0(J=2-1) lines (115.271 and 230.538 GHz, respec- 
tively), where the full half-power beam widths (HPBWs) 
are 45" and 23", respectively. Observations performed 

^ The Swedish-ESO Submillimetre Telescope is operated 
jointly by ESO and the Swedish National Facility for Radio 
Astronomy, Chalmers University of Technology 
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during 1988-1994 were made with a cooled Schottky-diode 
mixer receiver, while the subsequent observations were 
done with SIS receivers. All receivers operated in the 
single-sideband mode, and the average system tempera- 
tures, corrected for the atmospheric and antenna ohmic 
losses, were 660 K and 390 K for the CO(J=1-0) obser- 
vations in the first and second epoch, respectively, and 
340 K for the C0(J=2-1) observations. The backends 
were acousto-optical spectrometers, with channel separa- 
tions of 0.69 MHz (1.8 and 0.9 kms^^ at the two rest fre- 
quencies) and total bandwidths varying between 0.5 and 
1 GHz. 

We used the dual beam switching mode with a throw 
of Ril2' in azimuth. In dual beam switching the source is 
placed first in the signal beam and then in the reference 
beam. The two spectra produced are then subtracted to 
generate a spectrum with a minimum of baseline varia- 
tion. Great care was taken to ensure the quality of the 
data. No observations were done below the elevation of 
30° or above 84° (M83 passes within 1° from zenith at 
SEST), and pointing was regularly checked and updated 
using the SiO(t;=l, J = 2 — 1) maser line of the nearby 
AGB-star WHya (w3° from M83). Total pointing errors 
were typically less than 3", and data with pointing offsets 
suspected to be larger than 4" were disregarded in the sub- 
sequent data reduction (when possible, map points with 
such data were reobserved). The intensity scale was cali- 
brated using the conventional chopper wheel method, and 
the internal calibration errors in the corrected antenna 
temperatures, TJ, is within ±10% according to the SEST 
manual. In this paper we use the main beam brightness 
temperature scale Tmb, which is defined by Tmb = T^/ilmh- 
According to the SEST manual the main beam efficiencies 
[77mb(115 GHz)=0.7, ?7mb(230 GHz)=0.5] have been con- 
stant throughout the period of our observations. 

The observations were centered on the coordinates 
13'^36'"59!4 -29°52'05" (J2000), which is the optical cen- 
ter given bv lde Vaucouleurs et al ](|l97^.ForCO(J=l-0) 
we have 3634 spectra in 1900 positions with 11" spacing, 
and the coverage is complete out to a radius of 4'20". The 
grid extends parallel to the equatorial coordinate system, 
and the dense spacing (1/4 of the HPBW) was chosen 
in order to facilitate the use of deconvolution techniques. 
The integration time per position was chosen to render an 
rms around 70 mK (Tmb-scale) at the original velocity res- 
olution (1.8 km s^^). During normal weather conditions 
we had a typical on-source integration time of 60 seconds. 
Among the final 1900 spectra we have > 3a detections in 
1201 positions, i.e. a detection rate of 63%. The average 
rms noise in our data set is 74 mK, which should be com- 
pared to the peak temperature which usually lies around 
0.3 K (Tnib-scale) in the spiral arms in the outskirts of the 
map. 

In C0( J=2~l) we have two maps with different spac- 
ings covering two, partly overlapping, regions. The inner 
5'x3' is covered with a 7" grid, while the rest of the op- 
tical disk (out to a radius of 4.'5) is covered with an 11" 
grid. In total, we have spectra in 2574 positions, out of 



which we have > 3a detections in 1898 positions (73% de- 
tection rate). The spectra have an average rms of about 
90 mK at the original velocity resolution (0.9 km s~^). The 
observational parameters are summarized in Table |21 

3.2. Data reduction 

The data reduction was carried out with CLASS^, drp'^, 
and xs*. During the CO(J=1~0) observations, the num- 
ber of channels and the frequency resolution varied be- 
tween the spectrometers. We therefore resampled all data 
to 1.8 kms^^ (0.69 MHz) and by selecting the central 700 
channels this gave us a (fully sufficient) total bandwidth 
of 1260 km s-^ 

After flagging bad channels, and replacing them with a 
value obtained from interpolation between adjacent chan- 
nels, all spectra at the same position were averaged using a 
weight determined by the noise level. The baselines were 
stable, and only linear polynomial fits were subtracted. 
From the spectra we created a fits data cube for each 
transition, where the first two axes represent spatial coor- 
dinates, and the third axis represents frequency (or veloc- 

ity). 

We adopted a relatively dense grid spacing, compared 
to the beamwidth, in order to be able to increase the 
spatial resolution of the data using a MEM-deconvolution 
routine (maximum entr opy method) ; th e Statistical Image 
Analysis (SIA) routine l)Rvdbeckll2000|) . This procedure is 
particularly effective on a data set as densely spaced as 
ours (spacing in the range 1/4-1/2 of the HPBW), due 
to the high redundancy. The routine starts by deconvolv- 
ing a velocity-integrated intensity map, where the velocity 
range is chosen to be large enough to encompass all the 
emission. Once this is done, the procedure splits the initial 
velocity range into three equally wide regions, calculates 
their velocity-integrated intensities, and then deconvolves 
these three maps, using the deconvolved map from the 
previous step as the input guess. The deconvolution con- 
tinues in this hierarchical manner, successively narrowing 
the velocity interval, until a velocity range of 5kms~^ is 
reached. To deconvolve beyond this point is not mean- 
ingful since the S/N-ratios in the (narrow velocity range) 
maps drop below a critical value. For each of the two CO 
transitions we created 81 (3^) maps per transition (each 
covering about 5kms~^) which we assembled into a data 
cube per transition. By convolving the MEM-deconvolved 
data cube by a Gaussian with a FWHM equal to the ob- 
servational HPBW and comparing with the raw spectra, 
we verified that the MEM-results are reliable, and that the 
raw data set is homogeneous, i.e., free from variations in 

^ Continuum and Line Analysis Single-dish Software from 
the Observatoire de Grenoble and Institut de Radio- 
Astronomie Millimetrique 

^ Data Reduction Package developed at Onsala Space 
Observatory 

^ XS: Spectral Line Reduction Package for Astronomy devel- 
oped at Onsala Space Observatory 
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CO(J=1-0) 


C0(J=2-1) 


Center of the CO map 


R.A. = f3''36™59^4 (J2000) 




Decl.= — 29°52'05' 


' (■J2000~) 


Rest frequency 


115.271204 GHz 


230.53799 GHz 


Beam size 


45" 


23" 


Center velocity (LSR) 


520 km s"^ 


520 km s"^ 


Average rms (Tmb scale) 


74 mK 


90 mK 


Channel separation 


<0.69 MHz 


0.69 MHz 


Number of channels 


>700 


1000 


Main beam efficiency (77mb) 


0.7 


0.5 


Aperture efficiency (t^a) 


0.58 (27 Jy/K) 


0.38 (41 Jy/K) 



gain and pointing errors. The angular resolutions of the 
MEM-cubes are estimated to be «22" and «14" for the 
CO(J=1-0) and C0(J=2-1) data, respectively, and the 
velocity resolution is 5kms~^. 

We have also constructed convolved data sets by con- 
volving the raw data cubes with a HPBW beam size of 
20" and 15" for the CO(J=1-0) and C0(J=2-1) data, 
respectively (the spacing in these cubes are 11"). Even 
if this degrades the resolution by wlO-20%, the average 
rms noise in these cubes decreases from 74 mK to 25 mK 
and from 90 mK to 31 mK, respectively. The dramatically 
decreased noise level is due to the fact that the spacing 
is small with respect to the HPBW of the convolution 
kernels. Since the peak intensity in the spectra are rel- 
atively unaffected by this process, the average S/N-ratio 
in the spectra increases by a factor of w3. As a result of 
this, the detection rate increases to 92% and 89% for the 
CO(J=1-0) and C0(J=2~1) data, respectively (cf Sect. 
I3.1|l . These cubes were used to cross-check the results we 
got from the MEM-deconvolved data. Some results were 
also derived directly from the convolved cubes, especially 
in cases where the velocity resolution was more important 
than the spatial resolution. Furthermore, we convolved the 
raw C0(J=2-1) cube with a beam size of 44" and regrid- 
ded it to 11" spacing in order to compare with the con- 
volved CO(J=1-0) cube at the same resolution (Ri49"). 

In order to avoid problems with baseline irregularities 
when calculating the zero moment (over velocity) of the 
intensities we used a sliding window technique where at 
each position integration was performed inside a velocity 
window whose center and width depend on the position 
within the galaxy. We checked carefully that the window 
always covered all the emission. 

3.3. Correction for the error beam 

A closer inspection of the C0(J=2-1) spectra reveals, 
apart from the expected main component seen in the 
CO(J=1-0) spectra, a broad (velocity width 400 km s^^) 
and faint (peak intensity around 0.03 K) component. 
Despite the low peak intensity of this feature, its in- 
tegrated intensity is quite appreciable due to its large 
width. Typically, it makes up to 20%~35% of the total 
integrated intensity. This feature was discovered in the 



MEM-deconvolved data, which are noise-free, but it can 
be seen also in spectra produced by convolving the raw 
C0( J=2-l) data cube with a 45" beam. The appearance 
of this broad component changed with position in the disk, 
and we found that its shape could, in any given position, 
be reproduced very well by convolving the data set with 
a relatively large Gaussian beam. This indicates that the 
SEST picks up emission from a much larger area than the 
HPBW of 23", i.e., at 230 GHz the SEST suffers from a 
non-negligible error beam. 

We have corrected the C0(J=2-1) data sets for this 
error using a method described in Appendix We 
have found no trace of an error beam contribution in 
the CO(J=1-0) data (not unexpectedly, considering the 
higher main beam efficiency), and consequently we have 
not applied any correction to these data. 



4. The data 

4.1. Spectra 

In general, the spectra have a smooth and symmetric 
shape, but in some areas they show deviations. The spec- 
tra close to the nucleus are asymmetric, show extended 
wings, and a few of them have multiple peaks. The main 
reason for this is a large velocity gradient over spatial 
scales small compared to the telescope beam. Further out 
in the disk there are a few (interarm) regions which show 
multiple-peak spectra, such as the "Goul d-belt structure " 
located about 3' SE of the galaxy center l)Comer6nll207il]) . 
However, in the vast majority of the observed positions the 
spectral feature can be well described by a single, Gaussian 
profile. 

4.1.1. The CO(J=1-0) spectra 

In the arms we find that the peak temperatures are typ- 
ically 0.2-0.4 K, and in the inter-arm regions they are 
about O.IK. The spectra are widest in the nuclear region 
with (Gaussian) widths of 80-100 km s~^. Further out the 
linewidths are of the order of 7-14 km s^^ in the arms, 
with no systematic difference between arm and inter-arm 
regions. 
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The global CO(J=1-0) spectrum (left panel in Fig.Pl 
has a velocity-integrated intensity of 8.18±0.04Kkms~^ 
(excluding the calibration error) and a noise level of 
1.8 mK (all errors in this paper are 3a unless other- 
wise noted). The receding part of the galaxy (SW) has 
more CO(J=1-0) emission than the approaching part 
(NE): the velocity interval 510-610kms-i is 18% more 
luminous than the 410-510kms~^ interval. An excess of 
high er- velocity emission has also been se en in the atomic 
gas. iHuchtmeier fc Bohnenstengell l|l98lj) found that the 
higher-velocity part of the HI emission is 16% more lumi- 
nous than the lower- velocity part. 

The peak intensity spectrum (Fig. ^ right panel) has 
been obtained by selecting, for every channel, the highest 
intensity in the map. Therefore, the noise level in this 
spectrum reflects the noise in the overall worst case. From 
studies of the channel maps we conclude that the peaks 
in this spectrum arise from the nuclear region (487, 514, 
and 534kms~^), the western bar end (580kms^^), the 
eastern bar end (460 km s^^), and the NE outer spiral arm 
(437kms-i). 

By adding the CO(J=1-0) spectra (from the raw data 
cube) in concentric, 20" wide annuli (compensated for in- 
clination and position angle and centered on the IR nu- 
cleus, see Table^, we produced the spectra in Fig.|31 They 
show how the average C0( J=l-0) velocity-integrated in- 
tensity changes with galactocentric radius (which is indi- 
cated, in arc-seconds, in the upper right corner of each 
panel). The coverage in the map is complete to a radius 
of 229". 



4.1.2. The C0(J=2-1) spectra 

The global and peak intensity C0(J=2~1) spectra show 
characteristics very similar to those of the CO(J=1~0) 
data. Fig. 121 [since the C0(J=2-1) spectra were obtained 
with two different spacings, we regridded these data to 
a common 11" spacing before producing the spectra in 
Fig. 12]. The main difference is that the peak intensity 
spectrum is brighter than the corresponding CO(J=1-0) 
spectrum. This is most likely an effect of the emitting ob- 
jects being smaller than the CO(J=1-0) beam. 

In Table ^ we summarize the characteristics of the 
global spectra. The global line intensity ratio R21 = 
J T^b{2 - l)dv / J T^b{l - 0)dv is 0.77±0.01, based on 
the errors from Table 0] Note that, including errors in 
the calibration and the main beam and error beam ef- 
ficiencies, the result is 0.77±0.16. However, it should be 
noted that the relative calibration of the C0( J=l-0) and 
C0(J=2-1) data may be better than the absolute cali- 
bration and hence the error in their ratio can be smaller 
than indicated by the absolute error. 



Table 4. Peak intensities, velocity-integrated intensities 
(with errors), intensity- weighted velocities, and equivalent 
widths of the global CO(J=1-0) and C0(J=2-1) spectra 
in Figs. Hand 12 The errors are estimated from the noise 
levels. In the column for integrated intensities we also show 
the absolute errors. 



Peak Integral Centroid Eq. width 

[K] [K km s~^] [km s~^] [km s~^] 

CO(J=1-0) 0061 8.18±0.04±0.82 515±1 134±1 

CO(J=2-l) 0.046 6.28±0.03±0.62 513±1 135±1 



4.2. The maps 
4.2.1. CO(J=1-0) 

Figure 0] shows maps of the CO(J=1-0) velocity- 
integrated intensity (hereafter /i_o) and the peak inten- 
sity (hereafter Ti_o) for both the convolved data set (left 
column) and the MEM-deconvolved data set (right col- 
umn). The CO(J=1-0) emission is concentrated to the 
nucleus, the central part of the bar, the bar ends, and the 
spiral arms. The velocity-integrated intensity at the nu- 
cleus and the eastern and western bar ends are: 73, 24, 
and 30Kkms~^, respe ctively, which is ve r y clos e to the 
corresponding values in ICrosthwaite et "all |2002|) : 73, 22, 
and 28Kkms~^. In the MEM-deconvolved data set these 
features are brighter due to the increased spatial resolu- 
tion: at the eastern bar end we measure 49 K km s~^ , at the 
western bar end 52Kkms~^, and at the nucleus the max- 
imum /i_o is 235Kkms-i (ri_o=2.8K). By fitting an 
ellipse to the MEM-deconvolved /i_o-distribution we find 
that the CO peak coincides with the IR nucleus within 
the absolute position accuracy of our maps (3"). The axes 
ratio of this ellipse is 60": 29" with a position angle of 
36°. In the MEM-deconvolved maps the spiral arms are 
resolved, and they can be followed for almost 360°. The 
emission in the spiral arms breaks up into regularly spaced 
maxima, which we interpret to be giant molecular associa- 
tions (GMAs), a p henomenon previously seen in e.g. M51 
ijNakai et al.lll994l) . The outer arm, to the SE of the nu- 
cleus, seems somewhat disturbed in the sense that in the 
MEM-deconvolved Ti_o map it is possible to discern two 
parallel arms. Further toward the west, this arm appear 
to almost merge with the inner arm, which also appears 
to be the case on the opposite side of the galaxy (NE of 
th e nucleus). 

iRand et alJ |]J99l) (hereafter RLH) mapped the 
CO(J=1-0) emission in the inner eastern arm using the 
OVRO Millimeter Array (spatial resolution 6" 5 x 3"5). In 
Fig. |31 we show their data overlayed on our deconvolved 
/i_o map. The components close to the nucleus (compo- 
nents 1-5 using their nomenclature), are not prominent 
in our data. They are distributed along the leading edge 
of the bar. The reason is probably that a shock in this 
region creates substructures (in temperature, in density, 
and/or in velocity) in the CO distribution which elsewhere 
is smooth and therefore undetectable by an interferometer 
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Fig. 2. Left: Global C0(J=2-1) spectrum. Right: Peak intensity C0(J=2-1) spectrum. Before producing these 
spectra the C0(J=2-1) data cube was regridded to 11" spacing. 



(RLH found that only 2-5% of the single-dish emission 
was recovered in their interferometer observations). Our 
eastern bar end is located between components 6 and 7, 
and the spiral arm segment which continues southward is 
located between components 9 and 10. 

4.2.2. C0(J=2-1) 

The C0(J=2-1) maps show structures very similar to 
those in the CO(J=1-0) maps, however with a higher an- 
gular resolution, Fig.|S| The deconvolved map shows 
the "double nucleus" previously seen in the C0( J = 3 — 2 ) 
and CO(J = 4-3) lines by Petitna,s fc WilsonI 
This "double nucleus" is apparent also in our raw 
C0( J=2-l) data, but it is not possible to see in the con- 
volved /2-1 map, since this map is made using data which 
was both convolved (with a 15" beam) and regridded (to 
11" spacing). Since the gray-scale and contour levels in 
Figs. 0] and El are the same, we can directly infer that the 



C0( J=2-l) emission is less concentrated to the arms than 
the CO(J=1-0) emission. 

5. Results and discussion 

5.1. Molecular mass in M 83 

The H2 column density in external galaxies is commonly 
derived from the velocity-integrated C0( J=l-0) intensity 
(/co = / Tmhdv) using the expression 

7V(H2) = Xco X /co , (1) 

where Xqq is a conversion factor. The actual value of 
this factor, and its applicability in different environ- 

I '. 1 ' K 

ments have been wide l y discussed, see e.g. Wilson (1995); 
lArimoto et all l)l996() : ICombesI tOOd) : lYoung fc Rownd 
l|200l|) . In this paper we ha ve adopted the val ue Xco= 
2.3x10^° (Kkms-i)-i cm-2 ||Strong et al.ll988|) . which is 
supported by the analysis done in RLH. In order to con- 
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Fig. 3. C0( J=l-0) spectra produced by averaging all spectra within successive, concentric (and inclination-corrected) 
annuli. The average radius (in arc-seconds) of each annuli is shown in the upper right corner. The center spectrum is 
shown in the upper left corner, and the outermost spectrum in the lower right corner. Note that the intensity scale 
varies from panel to panel (the scaling factor is shown in the upper right corner of each panel). The velocity resolution 
has in this image been reduced to 5 km s~^. 



vert the measured intensities into mass surface density we 
used the expression 



I](H2) = 3.7 X Ico X cos z = 3.3 x Ico Mq pc" 



(2) 



where i is the inclination of the galaxy, and Iqo is given 
in Kkms~^ (This expression can be obtained by integrat- 
ing Eq. Using this conversion factor to convert the 
integrated intensity of the global CO(J=1-0) spectrum 
(8.18Kkms~^, Tab. 01 we find that the average (over 
the entire optical disk, and including the bar and the 
nucleus) surface mass density is 27.2 pc~^, and that 
the total molecular gas mass within the surveyed area is 
3.3x10^ Mq. 

5.2. The radial distribution 
5.2.1. Molecular gas 

We have, using the sliding-window technique described 
above, calculated the velocity-integrated intensity, cen- 
troid velocity, peak intensity, and center and width of 



a Gaussian emission profile fitted to the spectrum for 
all spectra in all data cubes. The left panel of Fig. [T] 
shows the C0( J=l-0) velocity-integrated intensities (ob- 
tained from the convolved data set) as a function of the 
(inclination-corrected) distance from the galaxy center. 
The scatter is dominated by the arm-interarm contrast. 
The data in the right panel of Fig. |7| were obtained from 
the integrated intensities of spectra averaged within con- 
centric annuli (see Fig. O and Table by converting the 
CO(J=1-0) integrated intensity into mass surface density 
as described above. In this procedure we used the raw data 
set, since we needed neither the higher spatial resolution 
of the MEM-deconvolved spectra nor the lower noise of the 
convolved spectra. The error bars reflect the uncertainty 
in the azimuthally-averaged radial distribution. To these 
data we fitted an H2 mass surface density distribution 
which is a combination of a Gaussian component and an 
exponential function (for the disk). 
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Fig. 4. The CO(J=1-0) maps. Velocity-integrated intensity in the convolved data set (resolution Ri49"; upper left) 
and in the MEM-deconvolved data set (resolution «22"; upper right). Contours are: 1, 3, 5, 7, 9, 13, 18, 27 (first white 
contour), 37, 50, 100 and 175Kkms~^. Maximum intensities in the convolved and MEM-deconvolved maps are 127 and 
235Kkms~^, respectively. Peak intensities in the convolved data set (lower left) and in the MEM-deconvolved data set 
(lower right). Contours are: 0.05, 0.1, 0.2, 0.3, 0.5 (first white contour), 0.7, 0.9, 1.2, 1.5, 2.0 and 2.5 K. Peak intensity 
at the center are 1.4 and 2.8 K for the convolved and deconvolved set, respectively. 



120 exp 



2265 pc 



Mq pc- 



(3) 



The Gaussian shape of the inner part of the distribution 
depends mainly on the convolution of the true distribu- 
tion with the beam (the scale length 649 pc in Eq.O cor- 
responds to a FWHM of 50"). The right panel of Fig. □ 



shows also the mass surface density of the HI gas. These 
data were obtained from Tilan us & Allen ( 1993() (here- 
after TA). The column densities in TA were obtained with 
the VLA, and the authors noted that their estimated to- 
tal HI mass within the Holmberg radius (0.87 x 10^ M©) 
was lower than the corresponding masses estimated by 




other authors using single dish observations, e.g., both 
iRoasta d et all l)l974l) and iHiichtmeier fc Bohnenstengei 
1,1981.) obtained about 1.6 x 10^ Mq (all masses recal- 
culated to the distance 4.5 Mpc). The most probable ex- 
planation for the deficiency was, according to TA, the ab- 
sence of short-spacing information. In order to compensate 
for this we have added the "missing" mass surface density 
(2.4 Mq pc-2) to the HI data of TA. The resulting HI mass 
surface densities are still low and typically H2 gas consti- 
tute about 70% of the total mass of neutral hydrogen (H2 
-|- HI) in the disk inside 6.8 kpc. 

Based upon the H2 mass surface density distribution 
(Eq. 121) we estimate a total molecular gas mass within a 
radius of 7.3 kpc (the extent of the map) of 3.5 x 10^ M©, 
and extrapolated to the Holmberg radius (7.'3 or 9.6 kpc) 
the mas s becomes 3.9 x 10^ Mfr^ (w hich agrees well with the 
result of lCrosthwaite et all l)2002|) . 3.7x 10^ M© (cf. also to 
the result 3.3 x 10^ Mq that was derived in Sect. I5.ll using 
the velocity-integrated intensity in the global spectrum). 
This is more than twice the total HI mass within the 
same galactocentric radius. It should, however, be noted 
that 80% of the HI mass in M83 is found outside the 
Holmberg radius l|Huchtmeier fc Bohnenstengeilll98ll) . As 
a comparison, t he ste llar mas s of M 83 is ^ 37x 10^ Mq, 
given V§.=7.37 l|de V aucoulcu rs et al.ll99^ and assuming 
a mass-to- light ratio of 2. 

In the Gaussian/exponential disk decomposition the 
major part of the molecular gas mass within the Holmberg 



radius lies in the disk, 3.6 x 10^ Mq. iLordl l)l987|) made 
CO(J=1-0) observations in 21 positions in the central 
parts of M 83, and he estimated the total molecular gas 
mass within 115" to be 1.2x10^ Mq (we have compen- 
sated for the difference in Xco and distance). Our H2 
distribution gives 1.4x10^ Mq. 

5.2.2. Comparison with data at other wavelengths 

In Fig.|Hlwe compare the CO radial distributions with cor- 
responding photometry data in the B, K, and Ha b ands. 
The B and Ha data come from Talbot e t all l|l97ij^. and 
the K data from.Adamson ct al. ( 1987). The scale lengths 
of the various tracers compare well over most of the radial 
range, but the comparatively low spatial resolution affects 
the CO curves in the central region and at the bar ends. 
The latter are responsible for the peaks at '^120". We 
have fitted a single exponential to the data in the radial 
range 140"-240", which corresponds to the region outside 
the bar but inside the area where our CO survey is com- 
plete. The scale lengths are: 135±14", 126±16", 112±20", 
100±14" and 132±35" for CO(J=1-0), C0(J=2-1), K, 
B, and Ha, respectively (the errors of the latter three do 
not include the measurement errors). Thus, the disk scale 
lengths are, within the errors, the same for these star for- 
mation tracers. 

The CO intensities decrease exponentially over the en- 
tire surveyed area, and we do not observe the sharp cut-off 
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Fig. 6. Same as in Fig.^jbut for C0( J=2-l) . The contours are also the same. The resolution of the maps are 27" and 
12" for the convolved and the MEM-deconvolved data set, respectively. For the convolved data maps the values at the 
center are llOKkms"^ and 1.3K. In the MEM-deconvolved map, the nucleus breaks up into two discrete sources, with 
similar velocity-integrated and peak intensities: 240Kkms~^ and 3.0 K. 



at 4f5 seen bv ICrosthwaite et al.l l)2002l) . In fact, we find than in the intcrarm regions. Assuming a direct propor- 
CO emission with Ji_o ~ 1 Kkms~^ in selected spectra as tionality between CO emission and mass surface density, 
far from the nucleus as 5 '5. this agrees well with the lower limit of the arm~interarm 

mass surface density ratio, 2.3, found by RLH using purely 

kinematical arguments. 

5.3. Brightness distributions 

In Fig. 121 we show the MEM-deconvolved CO(J=1-0) 
We find that in the MEM-deconvolved map the average velocity-integrated intensity map superposed on a true- 
CO(J=1-0) intensity in the arms is about 2.5 times higher color image. This image is a composite of three images 
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Fig. 7. The C0( J=l-0) velocity-integrated intensity versus galactocentric distance is shown in the left panel (obtained 
from the convolved data set). The dispersion is dominated by the arm-interarm contrast, and the bar-ends are evident 
at ^2 kpc. In the right panel the H2 mass surface densities are shown with formal errors (the data are taken from 
the spectra in Fig. |3Jl, together with HI mass surface densities. Also shown, as a solid line, is the fit of an II2 mass 
distribution described in the text. It consists of a Gaussian central part (dashed line) and an exponential disk (dash- 
dotted line). 

Table 5. CO(J=1-0) and C0(J=2-1) velocity-integrated intensities of spectra averaged within inclination-corrected 
annuli. The annuli are fully sampled out to 230". The C0(J=2-1) data were regridded to the same spacing as the 
CO(J=1-0) so that the disk was sampled in the same way for both transitions. The difference in angular resolution in 
the data sets has an effect on the response at the center. The errors (Icr) does not include the calibration uncertainty. 
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in the filters B, V, and R, rendered as blue, green, and 
red, respectively. These images have been obtained at the 
ESO Danish 1.54 m t elescope by Soren Larsen and the 
data are described in iLarsen fc Richtled lll999l) . On the 
large scale the C0( J=l-0) emission follows the dust lanes, 
both in and between spiral arms, but at some locations 
the CO(J=1-0) emission is displac ed toward the region s 
where star formation takes place jWiklind et al.l Il990|) . 
Most notably, in the inner spiral arm on the eastern side, 
the C0( J=:l-0) ridge is located downstream of the dust 



lane with a separation of about 10-15" (also shown in 
RLH), as are the HI and H/3 emissions as shown by TA. We 
note that in general the OB associations straddle the rims 
of the CO(J=1-0) emission concentrations, presumably 
due to extinction in combination with the photodissocia- 
tion by the young stellar population. The outer and inner 
spiral arms appear to connect at several locations, and 
these connection regions often show an increased number 
of OB associations. The "Gould belt structure", at (1.'75,- 
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Fig. 9. The MEM-deconvolved velocity-integrated C0( J=l-0) intensity as contours superposed on an RGB-map pro- 
duced from images in B, V, and R. The optical images were obtained by Soren Larsen at the ESO Danish 1.5 m 
telescope on La Silla, Chile. Three stars have been used for image alignment. These are located in the lower left corner 
and in both upper corners and are marked with crosses. The inset shows the CO(,/=1-0) map in gray-scale. 



2.^25) in the coo rdinate system of Fig. is such a region 
l|Comer6nll2'00ll) . 

Figure E| shows the C0( J=2-l) velocity-integrated 
intensity as contours on an RGB image, where the colors 
represent star formation (Ha; red), infrared emission (I 
band; green), and dust lanes (filters V~I; blue). In the disk, 
the observed C0( J=2~l) brightness distribution strength- 
ens the conclusions drawn from the C0( J=l-0) emission. 
In the central region the higher angular resolution at the 
C0(J=2-1) frequency reveals additional details. There 
are two central nuclear components relatively symmetri- 



cally placed on opposite sides of the IR center. At these 
locations the dust lanes, on the leading edges of the bar 
attach to the nuclear ring seen in J-K bv lElmegreen et all 
(Il998l) . Their color difference map ( V-I) shows an area of 
high dust extinction at the location of the NE compo- 
nent, whereas at the SW component there appears to be 
relatively small amounts of dust. It appears as if the NE 
CO component lies above the disk (or the bulge) and the 
SW CO component li es below or behind the sarn e. This 
was also concluded bv lSofue fc Wakamatsul l)l994|) based 
on similar data. Along the bar, the C0(J=2-1) emission 
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Fig. 8. Azimuthally-averaged CO(J=1~0) and 
C0(J=2-1) velocity-integrated intensity (cf. Table 
in annular zones compared to surface brightnesses in 
the Ha, B, and K filters. Note how the bar ends appear 
as a bump at 120" in the distributions. The spatial 
resolution of the C0(J=2-1) is better than that of the 
C0( J=l-0) data, which results in the difference in slopes 
of the inner parts of these distributions. The dashed lines 
show the fits of exponentials discussed in the text. 



follows the dust lanes on the leading edges. In the disk the 
CO emission traces the Ha emission very well, although 
there are regions with a clear anti-correlation, i.e., sub- 
stantial CO emission but none, or very weak. Ha emission, 
such as the peak at about half an arcminute west of the 
nucleus. 

Figure shows the C0( J=1"0) velocity-integrated 
intensity on a gray-scale image of the HI column density 
(data kindly provided by R.P.J. Tilanus and R.J. Allen). 
The angular resolution of the HI map is slightly better 
than that of the CO map. The CO emission (which is 
assumed to be linearly proportional to the H2 mass surface 
density intensity) and the HI column density follow each 
other closely, with one outstanding exception: the nucleus, 
where very little HI is present. At some places in the disk 
the H2 and HI maxima are displaced from each other, but 
there is no apparent systematic trend that H2 is outside 
HI or vice verse. In general, the separations are of the 
order of 3-9", which is small compared to the resolution 
in these maps (23"). 

Local minima in the HI and H2 gas surface densities 
also correlate well. At (50", 180") in the coordinate system 
of Fig. 111! marked minima can be seen in both H2 and 
HI. The total gas surface density in this area is about 3- 
6Mqpc~^, of which only wlMopc"^ is in the form of 
molecular gas. Such low mass surface densities should not 
give rise to any massive star formation, and indeed, the 
Ha image shows very little emission here. 

Figure ^1 compares the H2 (as probed by the 
C0(J=2-1) emission) and the HI gas in more detail, the 
resolution being fvlO". In the main disk the C0(J=2-1) 



emission and HI emission follow each other closely and 
several of the concentrations in the C0(J=2~1) emission 
are also seen in the HI distribution. Again the exception is 
I the nuclear region. The mass surface density of HI at the 
' + two nuclear components differ markedly: the SW compo- 
5 nent has a peak mass surface density of 6 Mq pc~^ , while 
1 at the center of the NE component HI is actually seen 
^ I in absorption against a nuclear continuum source, which 
f again speaks in favor of the idea that the SW and NE com- 
I ponents lie on different sides of the disk (or a bulge). The 
2 " mass surface density of HI at the center is very much lower 
I than that of H2, where the highest value in the /i_o map 
corresponds to wTSOMqpc"^. In the /i_o map it is not 
, possible to distinguish the separate nuclei due to the rela- 
tively low spatial resolution and the coarse grid. However, 
using the average 2-1/1-0 hne ratio (see Sect. I5.4|l . the 
mass surface density is ^llOOMopc"^ for each of these 
components. As the physical properties of the molecular 
gas in the center may differ from that in the disk, appli- 
cation of the standard conversion factor may lead to an 
erroneous result. Indeed, it is likely that at least part of 
the emission in the center comes from a diffuse and non- 
virialized component. This would lead to an overestima te 
of the molecular gas mass (jPapadopoulos fc Allen|[2000|) . 

The explanation for the nice correlation may be that 
the HI is mostly a dissociation product, as TA discuss. 
A smooth distribution of 'primordial' HI gas would not 
be detected in the interferometer observations (in which 
roughly 50% of the HI gas was detected) while HI pro- 
duced as an effect of dissociation of molecular gas in the 
neighborhood of star formation would be localized, and 
therefore detected by the interferometer. 



5.4. The CO (J =2-1) /(J =1-0) line intensity ratio 

We have created a map of the C0( J=2-l)/(J=l-0) line 
intensity ratio, -R21, at 49" resolution. The mean value 
over this map is 0.72±0.19, which is within the errors 
the value d erived from the velocity-i ntegrated intensities 
in Table ^ ICrosthwaite et all l|2002t) estimate a ratio of 
1.1±0.2, which we suspect is higher than ours because they 
have not corrected their C0(J=2-1) data for the error 
beam contribution (considering the main beam efficiency 
of the NRAO 12 m telescope at this frequency, wO.56, 
the error beam is likely to be comparable to that of the 
SEST). In the nuclear region and along the bar the ratio is 
0.83±0.04. The i?2i map shows also that on the arms, the 
line intensity ratio is lower than in the interarm regions. 
This is illustrated in Figure [TSl where i?2i is plotted as a 
function of radius in a 45° wide sector of the galaxy (posi- 
tion angle 90°-135°). In order to quantify and compare the 
ratio on the arms to that in the interarm regions, we have 
selected two sectors where the arms are well separated. 
The position angle of these sectors are: 90°-135° (SE) and 
270°-315° (NW). In these sectors we selected spectra be- 
longing to the inner, outer and inter-arm regions, and av- 
eraged them separately. Finally, "arm" -spectra were ere- 



A. A. Lundgren et al.: Molecular gas in the galaxy M83 



15 





Fig. 10. The MEM-deconvolved C0( J=2-l) velocity-integrated intensity as contours on an RGB image where the 
colors red, green and blue represent star formation (Ha), infrared emission (I band), and dust lanes (filters V-I), 
respectively. The inset shows the C0(J=2-1) map in gray-scale. 



ated by taking the average of the inner and outer arm 
spectra in the two respective sectors. From these eight 
spectra, four ratios were calculated. Assuming our error 
beam is correct (Sect. 0, the line ratio is significantly 
lower on the arms: 0.59±0.04 (SE) and 0.69±0.04 (NW) 
compared to the ratio in the interarm region 0.88±0.17 
(SE) and 0.98±0.21 (NW) (the errors are 3a and are es- 
timated from the noise level in the spectra). We believe 
that this is an effect of a (partly) different cloud compo- 
nent in the interarm regions, where the heating of the gas 
is more efficient than in the dense arm cloud complexes. 
As a note: We have also calculated the i?2i value for the 



MEM-deconvolved CO(J=1-0) and C0(J=2-1) maps at 
the common resolution 30". The arm value is roughly the 
same (0.6±0.1), but the interarm value is higher, 1.3±0.3. 

By averaging i?2i in concentric annuli, we find a trend 
in the line intensity ratio in the sense that it decreases with 
galactocentric distance. Fig. ^1 (this is also indicated in 
Fig.lHJ. In particular, the lower envelope of the i?2i-values 
is clearly decreasing. A similar difference between the R21 
line ratio in the disk and towards the nucl eus has previ- 
ously been shown to exist in our Galaxy ijSawada et alJ 
l200lh . The wiggles in Fig. El can be associated with fea- 
tures in the R21 map. The nucleus appears to be sur- 
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RA offset [^^J 

Fig. 11. The MEM-deconvolvcd velocity-integrated CO(J=:1-0) intensity as contours on an HI column density map 
(gray-scale) of roughly similar resolution («20"). The HI data comes from TA. The inset shows the CO(J=1-0) map 
in gray-scale. 



rounded by a ring of excited gas, which is reflected as a 
bump in the radial distribution. The bar ends and inner 
arms are reflected as a flat region at Ri 120-1 50" , while the 
low i?2i-values in the outer arms suppresses the average 
i?2i-ratio in the range 200-300". 



6. Conclusions 

We have mapped, in over-sampling mode, the C0( J=l-0) 
and C0(J=2-1) brightness distributions over the entire 
optical disk of M83. Significant corrections for an error 
beam have been applied to the C0( J=2-l) data. A MEM- 
method has been used to increase the angular resolution 
of the data. 



We find that: 
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RA offset [^^J 

Fig. 12. The MEM-deconvolvcd velocity-integrated C0(J=2-1) intensity as contours on an HI column density map 
(gray-scale) of roughly similar resolution («12"). The HI data comes from TA. The inset shows the C0(J=2-1) map 
in gray-scale. 



i) the CO emission is strongly peaked toward the 
nucleus, which splits up into two components in the 
C0(J=2-1) data. Also the bar ends are prominent. The 
CO emission follows the leading edges of the bar. 

a) Molecular gas spiral arms are clearly identified, and 
they trace, in most cases, the dust lanes. There are fre- 
quent "bridges" between spiral arms, and in some areas 
the arms are clearly disturbed. 



in) An average arm-interarm brightness contrast of 
about 2.5 is found for the CO(J=1-0) line. 

iv) The CO (J=2-l)/(J=l-0) line ratio, which is 
about 0.77 on average, differ between the arm and inter- 
arm regions. It is significantly higher in the latter. 

v) Regularly spaced molecular gas concentrations of 
mass ~1O''M0 lie along the arms. 

vi) The estimated total molecular gas mass is 
3.9xlO^M0, and within 7.3 kpc the H2 mass dominates 
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Fig. 13. The CO ( J=2-l)/( J=l-0) line intensity ratio, 
i?2i, as a function of galactocentric radius. The interarm 
region shows a clear trend of increased line ratio. The data 
are taken from a 45° wide sector SE of the nucleus (posi- 
tion angle 90°-135°). 




Fig. 14. The CO (J=2-l)/(J=l-0) line intensity ratio for 
individual positions (dots) and the azimuthally averaged 
ratio (with error bars) in annular zones in the plane of the 
galaxy. 



over that of HI by a factor of more than two. The esti- 
mated gas/stellar mass ratio is «0.1 in the optical disk. 

vii) The CO and HI emissions are very well correlated 
in the optical disk, and frequently concentrations of the 
two gas components coincide. 

via) The CO radial brightness distribution in the disk 
follows that of other starformation tracers as Ha emission, 
and continuum light in the B and K filters. The estimated 
scale length, of a fitted exponential, is about 120", corre- 
sponding to 2.6 kpc at the adopted distance. 
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Appendix A: Error beam 

The beam pattern of a radio telescope is affected by its 
surface accuracy. This means that the power received in 
the error beam may very well be comparable to the power 
received in the main beam. For sources that are compact 
relative to the extent of the beam, it is relatively easy to 
compensate for this underestimate of the true flux: The 
measured antenna temperature (T^) is simply divided by 
the main beam efficiency, T^mb, to attain the main beam 
brightness temperature (T^b)- The main beam efficiency 
is usually well known and can be found in the user man- 
ual of the different telescopes. For extended sources, the 
situation is more complex, since the antenna may receive 
power not only from the main beam. If the size of the error 
beam is comparable to the extent of the source, a simple 
multiplication of (?7mb)~^ will overestimate the flux at the 
position where the measurement is done. If the antenna 
pattern is well known, it is possible to calculate how the 
beam couples to the source and thus compensate for this 
effect. Unfortunately, the antenna pattern for SEST has 
never been measured properly, so this was not an option 
for us. Instead we devised a procedure to remove the im- 
pact of the error beam, using our own data set to estimate 
the strength and size of the error beam pattern. 

Nearly all C0(J=2-1) spectra did show two compo- 
nents: one narrow line, that agrees with the width and 
center velocity of the C0( J=l-0) spectra at the same po- 
sition and another, broader, component that is picked up 
by the error beam. By studying these broad wings in a 
number of spectra in the deconvolved C0( J=2-l) data set 
and by comparing with spectra in data cubes convolved 
with different FWHPs, we estimated that the size of the 
error beam at 230 GHz is of the order of 200". This is in 
good agree ment with the error beam size of 3', which was 
assumed bv I Johansson et all l)l998l|) . 

In order to try to estimate the impact of the error 
beam we produced a data set where we convolved the 
MEM-deconvolved data set with a 200" Gaussian profile. 
In each position an automated process, using the sliding 
window technique, calculated the velocity-integrated in- 
tensity in the region of the extended wings [i.e. outside 
the velocity range where we have seen emission in the 
CO(J=1-0) data] in both the deconvolved cube and in 
the 200" convolved cube. The ratio between these inten- 
sities corresponds to the fraction of the available emis- 
sion which is picked up by the error beam. We will refer 
to this as the error beam efficiency. In principle, the er- 
ror beam efficiency must lie in the interval 0.1-0.5, since 
the main beam efficiency is 0.5 and the moon beam effi- 
ciency is 0.9. Some values of our error beam efficiency esti- 
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Fig. A.l. These spectra show the original spectra (solid 
lines), 200" convolved spectra (short-dashed lines) and 
subtracted spectra (long-dashed lines) in four typical po- 
sitions. The coordinates with respect to the map center is, 
from upper left to lower right, (-42,-70), (-42,0), (-42,70) 
and (-42,140). The 200" convolved spectra are scaled with 
the derived error beam efficiency. 

mates are higher than 0.5. This happens in regions where 
the emission in the main beam falls outside the expected 
velocity range (e.g., in regions with exceptionally strong 
streaming motions or wide emission profiles), leading to 
an over-estimation of the flux in the wings. An attempt 
to make the window wider had a negative impact on the 
results since we were left with fewer channels to calculate 
the velocity-integrated emission in the wings. Neglecting 
values above 0.5, we found that the error beam efficiency 
was on average 0.27±0.11. 

Figure IXTl shows the original spectrum, the 200" con- 
volved spectrum (scaled with the error beam efficiency), 
and the resulting (cleaned) spectrum after subtraction of 
the two previous spectra, in four typical positions. The 
convolved spectra reproduces the shape of the wings very 
well, and the shape of the cleaned spectra resembles the 
corresponding CO(J=1-0) spectra (not shown). 

We have corrected the convolved and the MEM- 
deconvolved data sets for this error by subtracting the er- 
ror beam contribution for each individual spectrum. Based 
on the previous result, we took a conservative approach 
and fixed the error beam contribution to 0.27. In princi- 
ple, we could have let this ratio vary from spectrum to 
spectrum by fitting the 200"-convolved spectrum to the 
extended wings, but this would have introduced other, 
not easily controlled, errors. 
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